Cluster mass estimation from lens magnification 
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The mass of a cluster of galaxies can be estimated from its lens magnification, which 
can be determined from the variation in number counts of background galaxies. In 
order to derive the mass one needs to make assumptions for the lens shear, which 
is unknown from the variation in number counts alone. Furthermore, one needs 
to go beyond the weak lensing (linear) approximation as most of the observational 
data is concentrated in the central parts of clusters, where the lensing is strong. By 
studying the lensing properties of a complete catalogue of galaxy cluster models, 
one can find reasonable approximations about the lens shear as a function of the 
lens convergence. We show that using these approximations one can fairly well 
reconstruct the surface mass distribution from the magnification alone. 



1 Motivation 

A rich cluster of galaxies acts as a gravitational lens on the galaxy distribution 
beyond it. This simple fact can be used to derive a great deal about both the 
lensing cluster as well as the background galaxy population. Here we discuss 
how to best exploit the variation in galaxy number counts caused by the lensing 
cluster, which enables one to derive the lens magnification El. Recently, is has 
been shown that a depiction in number counts can clearly be observed for the 
clusters CI 0024+1654! and A1689§. 

However, in order to obtain a mass for the lens, or even a mass distribution, 
one needs the lens shear as well, because the lens magnification \x, shear 7 
and convergence k (the dimensionless surface mass density) are related as 

M" 1 = (l-«) 2 -7 a . 

There are ways for obtaining the shear from observations!], but one would 
like to obtain a mass estimate that is independent of other methods. This 
means we need to find assumptions for the shear, either as a function of surface 
mass density, or as a function oLthe lens magnification. We use a sample of 
numerical galaxy cluster models □ to find heuristic relations between 7 and k 
(or /x), some of which will have an underlying assumption about the physical 
state of the lens, like isotropy. 
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2 Estimating the lens convergence from lens magnification 



We wish to find a local relation for 7 vs. k and/or /1. The simplest one is the 
weak lensing approximation, fcu n = (1 — /i)/2, which is valid only in the linear, 
small-K regime El, i.e. in the outskirts of clusters. However, most observational 
data is available for the central parts of clusters, where this approximation is 
invalid, and we need to go beyond the weak lensing approximation to realisti- 
cally estimate the cluster surface mass density. 

There are only two local relations between 7 and k that result in, a single 
caustic solution of the magnification equation which is easily invertibleu: 7 = 0, 
corresponding to a sheet of matter, and 7 = K, for an isotropic lens. In the 
shearless case we have the estimate ko = 1 — 'PlfJ-l^ 1 ^ 2 , while for the isotropic 
case the estimator becomes Ki = (1 — 'P|/«| -1 )/2, where the V is the image 
parity, i.e. the sign of fi. Note that one can only measure and therefore V 
has to be assigned by hand. 

In. practice, substructure and asphericity of the cluster will induce extra 
shearLJ, especially in the surrounding 1ow-k neighbourhood, where substructure 
is relatively more dominant and filaments make the cluster most aspherical. 
An approximation that tries to take these cluster lens features into account, 
while still giving an invertible h(k) relation, is 7 = \/(c + c~ x — 2)k. This 
results in the amplification relation = (k^t c)(k — c _1 ), with caustics at 
k = c and k = cT 1 . The solution for k is then til 
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where we have introduced a second parity S which is the sign of k— (c+c 1 )/2. 
Note that the 7 = approximation is recovered by setting c — 1 . 



3 Testing the convergence estimators on cluster models 

Although the feasibility of actually obtaining fi from real data is an interesting 
topic for discussion, the issue here is how to proceed from the measured magni- 
fication to other properties of the lensing cluster, notably its mass distribution. 
We therefore assume here that one can reliably measure the lens magnification, 
and investigate how well the approximations allow us to reconstruct the lens 
convergence from just this magnification map. This should show us the best 
possible result each approximation can provide us, and reveals systematics. 

Using the thin lens approximation 13, we produce maps of the lens con- 
vergence, shear and magnification for the cluster models. We then use the 
magnification only (with full knowledge of its parity, though), to reconstruct 
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Figure 1: True versus estimated K-maps for a numerical model of a rich cluster. The top- left 
panel shows the true k, the next four panels show the estimates for the linear, 7 = 0, 7 = k, 
and 7 cx k 1 / 2 approximations, while the bottom-right panel shows the difference of the latter 
estimate compared to the true distribution (i.e. top-left panel minus bottom- left panel). 
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convergence maps using the various assumptions about the shear, and compare 
these to the true convergence. This is shown in Fig. 1 for the most massive 
model cluster from the catalogue, as this one has the largest range of possible 
values for the convergence k. 

The linear estimator performs very poorly, as it just follows the magnifi- 
cation, including the caustics. It is only doing well for small k, as expectedu. 
The 7 = assumption produces an overestimate for the convergence for all 
regions of the cluster. The k — 7 estimator underestimates the mass in the 
central regions of the cluster, and (slightly) overestimates for k < 0.2. The 
7 ex t/k estimator clearly performs best. 



4 Discussion 



The fact that these estimators can reproduce the surface mass density reason- 
able well for model clusters does not guarantee that they work on observational 
data. The best performing estimator has the clear disadvantage that two pari- 
ties need to be set. More generally, one has to deal with many intrinsic sources 
of error associated with observed magnification maps. Because the number of 
background galaxies that can be used to construct a magnification map will be 
finite, there will be shot noise. Therefore one needs to smooth the distribution 
of number counts, or average over annular bins. Furthermore, the backgrond 
galaxies are clustered and have different redshifts, which produces extija un- 
certainty. However, these problems can be dealt with to some extent Erfl, and 
the magnification method, employing our strong lensing estimators, seems a 
feasible way to directly measure the dark matter distribution of galaxy clusters. 
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